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Motivated by cosmological first-order phase transitions we examine the nucleation and evolution 
of vacuum bubbles in non-vacuum environments. Non-standard backgrounds can be relevant in 
the context of rapid tunneling processes on the landscape. Utilising complex time methods, we 
show that tunneling rates can be notably modified in the case of dynamical FRW backgrounds. We 
give a classification of the importance of the effect in terms of the relevant dynamical time scales. 
For both the bubble nucleation and evolution analysis we make use of the thin-wall approximation. 
From the classical bubble evolution on homogeneous matter backgrounds via the junction method, 
we find that the inflation of vacuum bubbles is very sensitive to the presence of ambient matter and 
quantify this statement. We also employ inhomogeneous matter models (LTB) and models that 
undergo a rapid phase transition (FRW) as a background and discuss in which cases potentially 
observable imprints on the bubble trajectory can remain. 
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1. Motivation 



Given the cosmic no hair conjecture, as well as the very long lifetimes of metastable states on 
the landscape, cosmological tunneling has been studied mainly between pure vacuum states. How- 
ever, in the context of the string landscape there are recently proposed scenarios where tunneling 
can be catalysed. When very rapid tunneling occurs a vacuum bubble can nucleate and find itself 
on a non-vacuum background that had not enough time to evolve to a de Sitter spacetime. Rapid 
tunneling occurs for instance in chain inflation [[[]] where many coupled fields are at work providing 
a series of catalysed tunneling processes through many minima on the landscape, or in the context 
of resonance or DBI tunneling [j2|, [3j]. 

Here, we will study the nucleation and the classical evolution of vacuum bubbles on non- 
vacuum backgrounds in the thin-wall approximation; for a more ample treatment see [Hp. A back- 
ground we utilise for both analyses is a flat FRW model that will serve as the setting for ambient 
phase transitions, like reheating. For the nucleation part we study the effects of a power law inflat- 
ing FRW background on the tunneling rates, using a complex time path formalism. For the classical 
part we probe the effects of homogeneous (FRW) and inhomogeneous (LTB) matter backgrounds 
on the evolution of the vacuum bubbles by employing Israel's junction method []|]. Note that inho- 
mogeneous backgrounds can also be relevant in the context of resonance tunneling, cf. [|], ^, ^]. 



2. Tunneling Rates 

The semiclassical calculation of tunneling rates between states of pure vacuum differing only 
in their values for A is a well known and commonly used result, first carried out by Coleman and 
De Luccia (CdL) [g]. CdL tunneling produces a bubble of new vacuum whose interior, and this is 
the main result of the CdL calculation, is also a de Sitter spacetime and thus can be used to model 
inflation. Hence, spherically symmetric bubbles of new vacuum expand into the old vacuum and 
so a first-order phase transition (PT) can be realised. 

However, the CdL mechanism and its variations rely on the participating states to be in a pure 
vacuum. Here we want to relax the de Sitter symmetry of the initial state and see how this can 
affect the tunneling rate. Of course, it is not at all clear how to do a calculation of tunneling rates 
on completely arbitrary backgrounds, and so we diverge only little from de Sitter spacetime in that 
we are taking background Friedmann universes in this section which will affect the tunneling rate 
through their time-dependent Hubble rate. 

Assuming spherical symmetry, a flat FRW background, and neglecting gravitational backre- 
action, the dynamics of the bubble becomes a 1 + 1 dimensional problem with the effective action 



dr/ 



-^■ea*{r\)T*{r\)-4lza a 3 (ri)r 2 (tj) y/ 1 - (d„r(ij)) 



(2.1) 



where tj denotes conformal time, r the coordinate radius of the shell, e the latent heat of the vacuum 
and a the bubble surface tension. The scale factor enters the equations of motion and introduces 
an explicit time dependence and this will make the tunneling rate time dependent as well. 

In JT0| ] a compendium of how one can treat tunneling processes in time dependent setups has 
been presented. Here we want to apply these techniques to the example of a power law inflating 
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background, a = {t\i/t\) +a with a(T]i) = 1, and thereby analyse the effects of a time dependent 
background on the tunneling rate. Note that for small deformations a this is just an inflationary 
slow -roll solution with slow roll parameter —d t H/H 2 a. For a power law inflation we get, after 
some transformations, the following equation of motion from the action 

This is an equation for T](r), which is an analytic (in general) complex function of real arguments 
r ranging between and oo. We can solve it only numerically, taking into account the boundary 
conditions for the canonical momentum p (t] ) = as well as d f r\ (0) = 0. We have carried out a 
parameter study with regard to the relevant model parameters a, T]o/t]i and e/a and found that the 
following three characteristic time scales can be used to describe the physical behaviour: the time 
scale tn, given by the inverse Hubble rate H 1 (tj), on which the background (scale factor) changes, 
the time scale tfj on which H itself changes, given by \d t H/H\~ (higher order derivatives of the 
expansion rate vanish in power law inflation), as well as the light crossing time for the bubble f cr0S s- 
(A) in case that the crossing time is the smallest scale in the setting, i.e. ? cross -C we find that 

the Minkowskian result Im^Mink = fr f° r the tunneling rate is a good approximation. Then, if 
the bubble crossing time scale is not much smaller than the inverse Hubble scale, ? cross > there 
are two possibilities. (B) if the rate of change of the Hubble time is still the largest time scale 
tH, ?cross ^ tfi we fi 11 ^ that our resu lt for the tunneling rate in de Sitter spacetime 

ImS d s = -j^ sinh 2 - In (l + (3Ho/e) 2 ) (2.3) 

in the quasistatic limit H = H (t]o) is a good approximation. However if finally (C) the crossing 
time cannot be regarded as small w.r.t. any other time scale, f cross > t^, then the tunneling process 
does become susceptible of the background dynamics. In case (C) the changing rate of expansion 
can have a significant effect on the tunneling rate as our numerical analysis confirms - in fact the 
tunneling rate in this case is enhanced w.r.t. the quasistatic approximation. Trying to exploit the de 
Sitter result ( |2.3[ ) in order to approximate the non-trivial behaviour of case (C) one should, instead 
of fixing H = //(rjo), average the expansion rate over one bubble crossing time range prior to 
tunneling and then insert the averaged Hubble rate into (2.3), see Fig. [I]. 

Another background also possibly of interest for chain inflation is a radiation dominated FRW 
universe with cosmological constant. Due to the cosmic no-hair theorem in its far future this 
universe will be vacuum dominated. However, its early phase where it is radiation dominated is 
characterised by a rapidly decreasing Hubble rate and so the dynamics may have consequences on 
the tunneling rate similar to the power law inflationary case. One difference is the existence of 
a particle horizon. The horizon scale is ~ t, as are the two background time scales. Because of 
causality the bubble crossing time is bounded by the horizon scale and thus it cannot become larger 
than any of the background time scales. As a consequence the corrections from dynamics should in 
general be small. Also, in the background spacetime there exists an initial singularity. A numerical 
analysis shows that the rates for nucleation of bubbles of about horizon scale become sensitive to 
the evolution of the background scale factor up to the vicinity of the initial singularity. Thus, the 
rates for such bubbles become sensitive also to the details of reheating and even earlier physics: the 
classical horizon problem also occurs in the context of vacuum tunneling. 
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Figure 1: Considering a FRW background with power law inflation numerical values for the tunneling 
rate as a function of the slow roll parameter d t H/H 2 are given (black circles). The values are normalised to 
the quasistatic approximation, that is the limit of the result for the tunneling on de Sitter background (2.3) 
when fixing H = H (rjo) with rjo denoting the nucleation time. Three typical time scales are notable here and 
their comparison can help to categorise the dynamics of the tunneling process. These are the bubble light 
crossing time f croS s, the inverse of the background Hubble rate tn as well as its respective rate of change ty. 
Only when the bubble crossing time is not negligible w.r.t. both of the other two time scales the dynamics of 
the background becomes significantly important for the tunneling rate and neither the Minkowskian nor the 
quasistatic limit are good approximations. One can improve the quasistatic estimate by employing a proper 
time averaged (dashed red line) or conformally time averaged (dotted green line) Hubble rate in (2.3). 

3. Classical Evolution 



In this section we want to study the (subsequent) classical evolution of a bubble on dynamical 



backgrounds, cf. [jllj. Like in the previous section we assume the bubble wall to be thin compared 
to its size. The spacetimes of the de Sitter bubble and the matter background will be matched along 
a common spherically symmetric hypersurface with induced metric findy'dyi = — dz + R 2 d£l 2 , 
stress-energy Sy = —ohij and the bubble surface tension a known from the preceding section. As 
the dynamical background spacetime we employ the spherically symmetric and radially inhomo- 



geneous Lemaitre-Tolman-Bondi (LTB) solution [12]: 



d^ 



2 (rd r a(t,r)+a(t,r)) 2 dr2 ? 2 

l+2£(r) V ' ' 



(3.1) 



The free function E(r) can be interpreted as either a measure of the local spatial curvature or as a 
measure of the total energy per unit mass of a shell at given radius r. Including dust and vacuum 
energy we get the following Einstein equations 

2 



d t a 
a 



2E 

a 2 r 2 



2M A 

+ 3 and 



a 3 r 3 



8 Tip 



2d r M 



(3.2) 



a 2 r 2 (rd r a + a) 

Here, M(r), the total mass within a given shell of radius r, is another free function and a third 
free function (the bang time) arises from the integration of the scale factor a. But there is also a 
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remaining gauge freedom in r and, when d,-M > 0, one can rescale r such that the function M(r) 
becomes M(r) = ^Ar^, where A is a constant. 

For our setup the junction conditions yield the following equations of motion in the LTB frame, 

-(l + 2E)rd t a + ^(l+2E) (1 + 27)7 (fd- r a + a) d t r 

o t r = — — r— ,d t o = p—= (3.3) 

(rd- r a + a) (2E - 2V) , 1 1+2E - {rd- r a + af {d t rf 



where r(t) denotes the bubble trajectory in LTB coordinates, Y = (rd t a) 2 — 2E + 2V with an effec 
tive potential V = 



\. There is a geometrical constraint on the 



surface tension due to the glueing procedure, 4na < y ^j- + A +~ A - _ 

We first consider the homogeneous limit of the LTB background, which can be obtained by 
imposing a homogeneous initial density distribution as well as a vanishing local curvature. The 
evolution of a comovingly nucleated {d,f = 0) bubble of new vacuum on such background substan- 
tially depends on the pressure balance and so on the relation of surface tension to latent heat and 



dust density. The initial force balance can be obtained by differentiating (3.3), 



12 i 1 M+-A_ „ 2p\ 

d 2 r = ~(-±- 2ko-^- \ . (3.4) 

The vital difference to nucleation on vacuum is that the dust obstructs the expansion of the bubble, 
eventually leading to a contraction. One can understand this from the competing forces: on the 
one hand the bubble tends to collapse under its surface tension and on the other hand the pressure 
induced by the latent heat of the vacuum (with A + > A > 0) works to expand the bubble. On a 
vacuum dominated background (p < p vac = A + /87i) the pressure force is able to outweigh surface 
tension. We find however that in a matter dominated background the presence of dust effectively 
makes the pressure support incapable of upholding the expansion of the bubble, see Fig. ^ (left). 

Next we will study bubble evolution on the radially inhomogeneous LTB background. Our 
aim is to probe whether the transition of the domain wall through background inhomogeneities can 
disturb the bubble trajectory. To implement inhomogeneities in the LTB model we can use the local 
curvature function E(r) as well as the initial density profile po(r) = p(ta,r). Recalling the results 
of the homogeneous limit, we see that the bubble would hardly be able to arrive at the inhomo- 
geneities because of the dust in the background. A numerical study indicates that even relaxing 
the condition of a comoving nucleation, i.e. considering a bubble with d t f(t ) > 0, does not save 
the bubble from collapsing whenever there is enough dust in the background. Nevertheless, there 
exist combinations of po(r) and E[r) for which the bubble can reach ambient curvature inhomo- 
geneities. For a radially growing initial dust profile po(r) small nucleated bubbles find themselves 
in a vacuum dominated region and are able to expand initially, cf. Fig. ||] (right), and can so reach 
the curvature inhomogeneities. We numerically traced the trajectories of these bubbles: for (poten- 
tially physical) interior observers there appears no difference w.r.t. the corresponding trajectories 
in the homogeneous limit, cf. Fig. ||. Note, however, that the amplitude of the curvature profile is 
additionally constrained by demanding the absence of shell crossing singularities. 

Another interesting dynamical background is a FRW model containing a perfect fluid that 
undergoes a rapid PT (e.g. reheating) while the vacuum bubble propagates through it. A numerical 
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Figure 2: Left: comovingly nucleated vacuum bubbles on a flat FRW background with dust and cosmo- 
logical constant (homogeneous limit of LTB). Shown are the possible initial bubble trajectories (as seen by 
comoving FRW observers) in the parameter space of the governing force balance quantities: the dust density 
p, surface tension cr and the latent heat of the vacuum £ vac = (A+ — A_) cf. (3.4). The expansion of 
de Sitter bubbles that is generic in the case of a pure vacuum background is here seriously hindered as long 
as the environment is matter dominated (p > p vac = A+/8?r). The forbidden area is due to a geometrical 
constraint from the junction method. Right: an LTB model with inhomogeneous initial dust distribution 
does allow expanding vacuum bubbles. Assuming, say, po(r) °< r 3 vacuum bubbles that have large nucle- 
ation scales are already in a matter dominated environment and are prevented from expanding while smaller 
nucleated bubbles find themselves in a yet vacuum dominated region and may initially expand. 

study confirms that the PT of the background does change the trajectory of the bubble notably and 
that, unlike in the previous case, this can also be seen from the interior, see Fig. |J 

4. Conclusions 

Rapid tunneling processes on the landscape may give rise to the nucleation of bubbles of new 
vacuum on a diversity of backgrounds. This kind of first-order PTs can be relevant for cosmological 
scenarios, like e.g. chain inflation. Therefore the nucleation as well as the evolution of de Sitter 
bubbles, which can be toy models of inflation, on any non-standard backgrounds is worth studying. 

Concerning the nucleation of vacuum bubbles we have analysed the possible effect on the tun- 
neling rate with an explicitly time dependent power law inflating FRW backround. We found that 
the rate can indeed be enhanced w.r.t. the quasistatic or Minkowskian approximations, if the cross- 
ing time scale of nucleated bubbles is larger than the time scales characterising the background 
dynamics. We have also seen that in a radiation dominated background the particle horizon con- 
strains the nucleation scale of bubbles and that bubbles nucleated with about horizon scale become 
responsive of details of the cosmology close to the Big Bang. 

In order to study the classical evolution of a de Sitter bubble on non-standard backgrounds 
we used Israel junction conditions to model bubble propagation in a matter environment. We saw 
already in the homogeneous limit that dust dominated backgrounds prevent bubbles from growing. 
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Figure 3: Numerical evolution of a vacuum bubble on an inhomogeneous matter background (LTB). Ex- 
panding bubbles are realised in an inhomogeneous initial dust distribution and are then traversing an addi- 
tional curvature inhomogeneity of the LTB model - the curvature profile is shown in the upper left panel. 
Comparing bubble trajectories on homogeneous (red) and inhomogeneous (blue) background, a significant 
difference is only seen in the exterior frame {upper right), while in the interior coordinates the transition oc- 
curs noteless {lower left). The sizeable effect of the transition in the evolution of the bubble surface tension 
{lower right) is to be considered as an artefact of the junction method rather than a physical effect. 



Nevertheless, for the inhomogeneous LTB background, initially expanding bubbles can be found. 
We numerically studied the effect of ambient inhomogeneity on such bubbles and found that pos- 
sible disturbances of the wall trajectory would not be seen by inside observers. We also analysed 
the effects of a rapid PT in a FRW background (e.g. reheating). Our analysis suggests that in this 
scenario perturbations of the bubble wall are present also for observers inside the bubble. 



Within the context of primordial bubble collisions it has been pointed out [13, 14] that similar 



disturbances of the bubble trajectory could in principle be observable in the CMB. 
References 

[1] K. Freese and D. Spolyar, JCAP 0507 (2005) 007, arXiv:hep-ph/0412145. 

[2] S. Sarangi, G. Shiu and B. Shlaer, Int. J. Mod. Phys. A 24 (2009) 741, arXiv:0708.4375 [hep-th]. 

[3] S. H. Henry Tye, arXiv:hep-th/061 1 148. 

[4] D. Simon, J. Adamek, A. Rakic and J. C. Niemeyer, JCAP 0911 (2009) 008, arXiv:0908.2757 [gr-qc]. 
[5] W. Israel, Nuovo Cim.B44S10 (1966)1, Erratum-ibid. B 48 (1967 NUCIA,B44, 1.1966) 463. 
[6] P. M. Saffin, A. Padilla and E. J. Copeland, JHEP 0809 (2008) 055, arXiv:0804.3801 [hep-th]. 
[7] E. J. Copeland, A. Padilla and P. M. Saffin, JHEP 0801 (2008) 066, arXiv:0709.0261 [hep-th]. 



7 



lst-order phase transitions beyond the standard inflationary scenario 



Aleksandar Rakic 





w = l/3- 



• w = — 1 



0.4 0.6 

v/A_/3f 




0.0 
0.0 









w = — 1 — > w = 1/3 



0.8 
t) 0.6 
0.4 




:-!-> IV =1/3 



0.4 0.6 

v/A_/3f 



Figure 4: Numerical survey of the propagation of a vacuum bubble in an environment (FRW) concurrently 
undergoing a reheating PT or vice versa. Three different initial values <Jq (red, blue, green) are indicated 
respectively. The upper row shows the resulting bubble trajectories in the FRW frame: while in the reheating 
case the expansion of the bubble is reversed by the passage through the ambient PT, in the w— l/3^w= — 1 
case the contraction of the bubble is stopped by the transition. Different to the LTB case notable effects 
remain to be seen also in the interior frame (middle row). Analogously to the LTB case a pronounced effect 
on the surface tension itself is found (lower row) which is a consequence of the jointly used junction method. 
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